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Abstract 

We report on the first mathematically rigorous proofs of a transition to a giant vortex state of a 
superfluid in rotating anharmonic traps. The analysis is carried out within two-dimensional Gross- 
Pitaevskii theory at large coupling constant and large rotational velocity and is based on precise 
asymptotic estimates on the ground state energy. An interesting aspect is a significant difference 
between 'soft' anharmonic traps (like a quartic plus quadratic trapping potential) and traps with a 
fixed boundary. In the former case vortices persist in the bulk until the width of the annulus becomes 
comparable to the size of the vortex cores. In the second case the transition already takes place 
in a parameter regime where the size of vortices is very small relative to the width of the annulus. 
Moreover, the density profiles in the annulus are different in the two cases. In both cases rotational 
symmetry of the density in a true ground state is broken, even though a symmetric variational ansatz 
gives an excellent approximation to the energy. 

1 Introduction 

A superfluid confined in a rotating anharmonic trap, where the rotation speed can in principle be ar- 
bitrarily large, undergoes several phase transitions as the speed increases. At first the fluid is vortex 
free [23 13] but then quantized vortices emerge, eventually forming a vortex lattice [5J [TH1 H31 [U 13 E] 
that may persist even when the speed is so large that the centrifugal force creates a 'hole' with strongly 
depleted density in the middle of the trap [THl [HI HI [3Z3 ■ Above a certain rotation speed a transition to 
a giant vortex state takes place. In this state the vortices disappear from the annulus where the bulk 
of the superfluid is concentrated while a macroscopic phase circulation remains. This phenomenon has 
been studied theoretically by variational and numerical methods in the past [T71 HH1 HH1 HOI 12H HH1 HI] 
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but mathematically rigorous proofs of the giant vortex transition have been obtained only very recently 
[3TI fTTl [T2l IT31 IT4l 132"] . An experimental realization of this transition appears to be still out of reach 
although anharmonic traps have been available already for some time [3] |MJ [3TJ [37] . In the following we 
report on the main findings of this analysis with emphasis on [131 114) . 

2 Setting the Stage 

2.1 The basic many-body Hamiltonian 

The quantum mechanical Hamiltonian for N spinless bosons with an external potential, V, and a pair 
interaction potential, v, in a rotating frame with angular velocity O rot is 

N 

j—X l<i<j<N 

Here Xj € Mr and Lj = — btj x Vj is the angular momentum of the jth particle. The Hamiltonian can 
alternatively be written in the 'magnetic' form 

N 

j=l l<i<j<N 

with the vector potential 

A(x) = fl rot x x = £l rot re,j 

where r denotes the distance from the rotation axis and e$ the unit vector in the angular direction. This 
way of writing the Hamiltonian corresponds to the splitting of the rotational effects into Coriolis and 
centrifugal forces. 

2.2 Harmonic vs. anharmonic traps 

If V is a harmonic oscillator potential in the direction _!_ to fi ro t, i-e., 

V(x) = itt trap r 2 + l/ll(z) 
then stability requires f2 ro t < f^ti-ap- Rapid rotation means here that 

^rot ^ ^trap 

from below. On the other hand, if V is anharmonic and increases faster than quadratically in the direction 
_L to r2 lot , e.g. V(x) ~ r s + V^(z) with s > 2, then J7 rot can in principle be as large as one pleases and 
rapid rotation means simply f2 rot —¥ oo. 

These two cases are quite different both physically and mathematically. The former leads to an 
effective many-body Hamiltonian in the lowest Landau level of the magnetic kinetic energy term in (fTJ) 
and bosonic analogues of the Fractional Quantum Hall Effect (see [3HJ [351 HI US])- In the case of rapid 
rotation in an anharmonic trap, as considered here, it is usually sufficient to employ Gross-Pitaevskii (GP) 
theory for an effective description. We remark, however, that a small anharmonic term appropriately 
tuned can also lead to interesting modifications of the Quantum Hall states of harmonic traps 33 . 
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2.3 The Gross-Pitaevskii limit theorem 



The following basic fact about the many-body Hamiltonian ([T]) for N —> oo with Na and f2 ro t fixed, where 
a is the scattering length of the (repulsive, short range) interaction potential v was proved in [27] : 

There is (possibly fractionated) Bose-Einstein condensation in the ground state as N — > oo, and the 
wave function of the condensate ( "superfluid order parameter" ) is a minimizer (in general not unique) of 
the GP energy functional 

£ 3 G d[*]= / {i|(iV + A)M/| 2 + (^-^ r 2 ot r 2 )l*| 2 +27r7Va|vI/|n (2) 

with J |W| 2 = 1. The Gross-Pitaevskii partial differential equation 

{ (iV + A) 2 + (V - \tf ot r 2 ) + 47r7Va|*| 2 }* = ^ 

with the chemical potential fi is the variational equation corresponding to this minimization problem. 

The rigorous proof of this theorem is far from simple, as can be seen from the fact that a Hartree 
variational ansatz for the Hamiltonian (|TJ) (that would anyhow only lead to an upper bound) is meaningless 
if the interaction potential has a hard core. Even for 'soft' potentials a naive computation would not lead 
to |(3J) with the scattering length as parameter, but rather the integral J v (that gives only the lowest Born 
approximation to the scattering length). For the mathematical background of this and related results 
[28j may be consulted. A limit theorem that holds uniformly the parameters f2 rot and Na as N — > oo, 
but is restricted to the leading order, was proved in [5]. 

2.4 2D Gross-Pitaevskii theory in anharmonic traps 

The GP minimization problem has two parameters, f2 ro t and Na. We shall be concerned with phenom- 
ena that occur in anharmonic traps in the asymptotic regime where both f2 rot and Na are large. For 
convenience introduce 

e = (2irNa)- 1/2 

which is small if Na is large. (In appropriate units e is the 'healing length'.) 

For traps that are sufficiently elongated along the rotational axis (z-direction) the properties of the 
condensate are to a good approximation independent of z and we may consider a 2E0 energy functional 

£^m= f (i|(iV + A)^| 2 + (y-IO? ot r 2 )|*| 2 + ^) 

with a trap potential of the form (for simplicity) 

V(r) = kr s (3) 

with s > 2, k > 0. Here fi rot can be arbitrary large. 

The limiting case s —> oo corresponds to a 'flat' trap with fixed boundary at r = 1. The effective 
potential is then simply — ^fl 2 ot r 2 and the integration is limited to the unit disc in R 2 . A word of caution: 
The limit s — >■ oo can not be interchanged with the limits e — y 0, f2 ro t — > oo as discussed in Section 3.6 
below. 

The analysis of the GP minimizer is guided by the following heuristics: 

• A vortex, i.e., a zero of the wave function 4 r (x) = |v?(x)| exp(i#(x)) with an accompanying nonzero 
winding number of the phase factor, reduces the kinetic energy because the associated current 
~ V0(x) compensates partly the velocity field generated by A(x) = J7 rot x x. 

1 A 2D description is, of course, also appropriate in thin traps where the motion along the z-axis is 'frozen' |36| . 
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• A vortex causes also a change in the density, however, (mass is moved from the vortex core to the 
bulk) and this increases the interaction energy that depends on the density at the potential location 
of the vortex. The energy balance decides whether or not a vortex is favorable, and if that is the 
case, the size of the vortex core. 

• A vortex is the more costly the higher the density. At sufficiently high rotational velocities the 
compression due to centrifugal forces creates a 'hole' and the density in the bulk increases until, at 
some point, vortices become too costly. 

2.5 Scaling of the energy functional 

The effective potential (kr s — ^Q 2 ot r 2 ) has a unique minimum at r = (Q^ ot / (sk)) l ^ s ~ 2 \ Taking this as 
a length unit we obtain the scaled energy functional 

£ GP M= / {h 2 \(iV + flxe^j\ 2 + n 2 W(xM 2 +e~ 2 \^} 

JR 2 

where x = |x|, tt ~ ^+ 2)/(;i ~ 2) , and 

W(x) = (±x° - \x 2 ) . 
The scaled potential has a minimum at x = 1, independent of CI. 

3 Analysis of the GP Minimizers 

3.1 Critical velocities 

The basic facts for traps of the form ((3]) with 2 < s < oo can be summarized as follows. 
As f2 increases there are three critical velocities: 

• Cl' cl ~ | loge| marking the appearance of the first vortexJl 

• ^c2 ~ £ marking the appearance of a 'hole' due to the centrifugal forces. 

• f2 C 3 ~ e~ 4 marking the transition to a 'giant vortex'. 

For the first transition we refer to [5J [2H HJ ISO]- For Cl c i fl <C f2 C 3 the ground state energy is 
well approximated by assuming a triangular vortex lattice in the bulki In the limit e — > the vorticity 
becomes uniformly distributed with density fi [14]. For ft > f2 C 3 the bulk is free of vortices but a 
macroscopic circulation around the origin remains [131 114] . 

3.2 The vortex lattice regime 

The ground state energy for fi c i «!!< C 3 can be computed exactly to subleading order [T4"] : 
Theorem 1 (Energy between fi C 2 and f2 C 3). //e _1 < ft £~ 4 as e — > 0, then 

E GP =E TF + !miog( £ 4 fi)|(l + (l)). 
2 Here W ~ e- 4 ^ a+2 ^n iot . This scaling is more convenient than Q ~ Q^+ 2)/(s_2) for fi < 1/e. 

3 The reason why a triangular arrangement with hexagonal unit cells is optimal amongst regular lattices can be made 
plausible by appealing to an electrostatic analogy and Newton's theorem |10| : Hexagonal cells are as close to being circular 
as possible and thus have smaller multipole moments and lower interaction energy than other cells. 
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Here E TF denotes the energy without the kinetic term. Below D,2 C a similar formula holds (with a 
different scaling, £~ 2 /( s + 2 ) as length unit): 

Theorem 2 (Energy between fi c i and f2 C 2). If | loge| <C fi' < as e — > 0, i/ien 

^gp' = ^tf' + i n ,| iog( e anO|(l + o(l)). 

3.3 Vortices reduce kinetic energy 

The potential term ~ fl 2 and the interaction term ~ e~ 2 become comparable when 

n ~ e _1 . 

This is the order of the second critical speed f2 C 2 above which the centrifugal force creates a 'hole'. In 
the sequel we shall focus on rotation speeds around and above f2 C 2 which means that f2 > e _1 . In this 
regime the kinetic energy term ^|(iV + A)^ 2 is formally also of order 1/e 2 if f2 ~ 1/e. Its contribution 
to the energy is, however, of lower order, namely ~ il|loge|, because a lattice of vortices emerges as 
e — >• and reduces the kinetic energy as remarked in 2.4. 

3.4 The giant vortex regime 

Consider a variational ansatz for the wave function of the for 

i/j(x) = g(x) exp(if2i?) 
with a real valued function g, normalized such that J g 2 — l.This gives 

£ GP W= I {h\Vg\ 2 + ln 2 { X -x^) 2 g 2 + n 2 {^-\x 2 )g 2 + s- 2 g i }^£^[g]. 
Jr 2 

The unique positive minimizer g gv of £ gv is rotationally symmetric and we denote the corresponding 
energy by . 

The following results are proved in [T31 [H] . 

Theorem 3 (Energy in the giant vortex regime). There is a constant < CIq < oo such that for 
O = Oq £~ 4 with Qq > Qq the ground state energy is 

E GP = £gv + (|l og£ |9/2)_ 

Theorem 4 (Absence of vortices in the bulk). There is a constant c > such that for f2 = f2g with 
Clo > flo and e sufficiently small the minimizer ip GF is free of zeros in the annulus 

A = {*: |l-a?| < cfi-^^logel 1 / 2 }. 

3.5 On the proof of the GV transition 

The main issue is a precise lower bound to the energy. We restrict £ gv to the annulus A, obtaining a 
positive minimizer go. Define u(x) on the annulus by writing 

V> GP (x) = g (a:)u(x)exp(ifitf). 
The function u contains all possible zeros of tp GP in the annulus. 

4 For simplicity of notation we assume that fl is an integer which is justified since Q — > oo. 
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The variational equation for g leads to the lower bound 



E GF >E%+£ A [u] 



with a functional of Ginzburg-Landau type with g$ as weight: 



£aM= I g 2 {±\Vu\ 2 -B-J(u)+e- 2 g 2 (l-\u\ 2 ) 2 } 



J A 




For this purpose one writes g 2 B = V ± F with V = (— d X2 ,d Xl ) and a potential function F. Integra- 
tion by parts and estimates of F (this is the key point!) give 



leading to the lower energy bound. 

A consequence of this bound, combined with the variational upper bound E^ < is an upper bound 
on the interaction term for large fto: 



J A 

Together with the upper bound to the kinetic energy and standard inequalities this implies that u 
must be close to 1, in particular free of zeros. 

3.6 Comparison with the 'flat' case 

The flat case, s = oo, that is treated in detail in [HJ[l2], differs from the case s < oo in several respects: 

• The GV transition takes place at ft ~ e~ 2 | loge| _1 rather than ft ~ e~ 4 . 

• The density profile in the GV regime is of 'Thomas-Fermi' type in the 'flat' case, but for s < oo it 
is gaussian around x — 1. 

• The 'last' vortices before the GV transition have size ~ e 3 / 2 that is much smaller than the thickness 
of the annulus ~ e|loge|. For s < oo the size of vortices, ~ e 2 and the size of the annulus, 
~ e 2 | log el 1 / 2 , are almost comparable. 

The techniques of proof in the two cases are also by necessity different: While vortex ball constructions 
and subsequent jacobian estimates (see [35]) for the potential function are applicable for the 'small' 
vortices in a 'flat' trap they are useless for s < oo and new ideas are required. 

3.7 Circulation and symmetry breaking 

At low rotation speeds below the onset of the second vortex the GP minimizer has rotationally sym- 
metric density, but a vortex lattice clearly breaks the symmetry. On the other hand, the giant vortex 
variational ansatz, that gives an excellent approximation to the energy and circulation for J7o > ^Oi is 
an eigenfunction of angular momentum. A true minimizer does not have this property, however: 

Theorem 5 (Circulation and rotational symmetry breaking). In the giant vortex regime Qq > CIq the 
circulation of any GP minimizer is 2tt ft + 0(ft o |loge| 9/4 ) ; but no minimizer is an eigenfunction of 
angular momentum. 

Theis result holds both for s < oo and s = oo [TT]-|14j. 
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4 Summary 



The study of the GP equation for dilute Bose gases in rotating, anharmonic reveals a surprising rich 
landscape, both from the mathematical and physical point of view. Detailed analysis can be carried out 
in an asymptotic regime where both the coupling constant and the rotational speed are large. Among 
the results found are: 

• Energy asymptotics corresponding to a distribution of vorticity in a lattice of vortices for f2 c i <C 

fi < n c3 . 

• Emergence of a 'hole' with strongly depleted density above a critical rotation speed Sl c2 ~ e _1 . 

• Transition to a 'giant vortex' state above f2 C 3 ~ e~ 4 where the vortex lattice disappears from the 
bulk and all vorticity resides in the 'hole', creating a macroscopic circulation in the bulk. 

• Breaking of rotational symmetry, also in the giant vortex regime. 
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